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Modified Bess el Functions (SADOI Only) 

Closed 

300 

0-10 inclusive 

At least 7 decimal places 

Approximately 190 m.sec. 

This routine calculates the modified Bess el Functions 

which are defined as: 



- co v+2m 

v ~ m=0 v+2m . P/, , , -, \ 
2 ml J,^v + m + l / ) 



K (z) = » .* [I (z) - I (z)] 
v x ' 2 ^.in. v n -v v 



These are solutions of the modified Bessel equation 
of order *v which is 



2 d 2 w , dw / 2 , . 2v _ 

dz. 



The accuracy of the routine depends upon the accuracy 

of the argument. For this reason, double precision argument 

may be used. In cases where the argument, scaled by a 

suitable scaling factor, has sufficient accuracy in single 

precision, a single precision argument may be used. 

The true value of the argument must be greater than .000001. 

Values smaller than this will give incorrect results . 

To use a double precision argument, where the argument is 

-M 
scaled by 2 , store the MSP in location OF and ISP in 

the accumulator and enter with 



P 



p+1 



50 MF 

50 pF 

^T7 



NOTE: 



- 2 



where M 70, but M< 39. 

If a single precision argument is to be used, the 
accumulator must be clear upon entry. One form for 
entry may be 

' 50 MF 



P ♦ 1 [ 26 ( )' 

Again the MSP must be in location OF. The routine will 
leave the following in the locations shown; 

in location 7F l/l6 e~ x L^x) 

» « 8F 1/16 e x K (x) 

» " 9 f 1/16 e" x I Q (x) 

" « 10F 1/16 xe x K-^x) for x < 2 

» « 10F 1/16 e x K^x) for xy 2 

Notice the two forms for K, (x) depending on the size of 
the argument. 

The aim of the double precision argument is to allow 
the use of a wide range of arguments, all of which may 
be scaled by the same scaling factor. This means, for 

example, that if a range of values for the argument is 

12 13 

from .00005 to 2 , they may all be scaled by 2 v . The 

routine rescales each value down to the minimum value 

needed for eacho For small arguments scaled by large 

factors, a double precision argument is needed to insure 

the accuracy of the routine. 

This routine contains the subroutines; (R l) Square Root, 
(S3) Natural Logarithm, (R 5) Reciprocal, and (S 4) 
Exponential. These may be used by other parts of the pro- 
gram and may be entered in the usual manner since they 
are referred to in this routine as symbolic addresses. 



- 3 - 

The method used to obtain these functions may be found 
in "Mathematical Tables and Aids to Computation" 1956* 
Vol. 10, pp* 162-163. 
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.OCATION 



10 
11 
12 

13 
Ik 

15 
16 



ORDER j 

ooK(yii) j 

kO IF 
K5 170L 
k2 108L 
10 20F 
42 188L 
kl 189L 
26 237L 
50 IF 
L5 F 
00 IF 

36 6l 
26 9L 
kO F 
F5 I89L 
40 189L 

LO 188L 
36 9L 

26 4l 

L5 F f 
kO 5F 
L5 188L 
LO I89L 
kO 6F 
L3 6F 
32 30L 
kO F 
F5 F 
36 29L 
40 F 
F5 F 
36 16L 
26 109L 

L5 16% <: 

LO 5F 



1 



NOTES 



Plant Link 



Store M 



Place argument 
in AQ 



Rescale arg ; 
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TTT] 



y"by new factor M' 



X • 2 N 

Store M' 
Test size of 

M 1 and transfer 

to proper portion 

of routine 



J 

a 



Test for X = 2 



LOCATION 



ORDER 
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17 
13 

19 
20 
21 
22 
23 
2k 

25 
26 

27 
28 

29 
50 
31 
32 
33 
3^ 



kO F 
L3 F 
32 19L 
26 22L 
50 158L 
L5 19L 
46 9^L 
L5 86L 
42 72L 
26 48L 
L5 190L 
LO 5F 
40 F 
L3 F 
36 26L 
L5 F 
36 ^7L 
26 109L 
L5 202L <H 
40 9F 
L5 179L 
40 7F 
22 138L 
50 193L 
L5 5F 
10 IF 
22 31L 
L5 5F 
10 2F 
kO 5F 
L5 5F 
00 IF 
50 77L 
50 33L 
26 (S3) 
10 6f 



r x = 2 
X > 2 



Test for 
X = 3*75 

X = 3.75 

X < 3.75 
X > 3.75 



Form X/k 



V Form 1/32 In x/2 



LOCATION 
35 

36 

37 

38 

39 

1*0 

in 

42 

44 

45 
46 

47 

49 
50 
51 
52 



ORDER 



75 187L 
00 2F 
1*0 176L 
U 5F 
50 69L 
50 37L 
26 (S4) 
40 F 
50 E 
7 J F 
40 F 
50 F 
7 J F 
40 177L 
50 95L 
50 42L 
26 (R5) 
10 3F 
S5 99L 
50 F 
40 178L 
L5 6lL 
42 72L 
26 53L 
L5 68L 
42 72L 
LI 5F 
50 48L 
26 (S4) 
40 F 
50 F 
7J F 
40 F 
50 F 
7 J F 
40 177L 
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► Form e 



-x 



Form e 



-x 



f Form l/8 e x 



1 



*■ Form e 



-x 



/ 



LOCATIOB 
53 

5* 

55 

56 

57 
58 
59 
60 
61 
62 

65 
66 

67 
68 

69 



L5 28L 
42 6ZL 
50 191L 

L5 5F 
66 190L 

35 F 
40 F 
L5 37L 
46 65L 
L5 191L 
40 199L 
L5 200L 
40 201L 
50 F 

7J F 
40 IF 
40 2F 
50 73L 
50 2F 
7 J F 
L4 199L 
40 199L 
F5 201L 
40 201L 

36 F 
50 IF 
7J 2F 
40 2F 
F5 62L 
42 62L 
26 62L 
50 85L 
L5 199L 
L4 192L 
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► Form X/3.75 



/ 



^ 



Form Summation 



of terms 



for. small 



arguments 



Form l/l6 e" x I n (x) 



ixxmxw 


ORDHR 
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■■■ TP;. . 


kO 3? 




1 




$0 181L 




7 for snail arg. 


?i 


50 177L 
7J-5P 






72 


40 JF 
26 F 


* 


/ 

1 


73 


L5 70L 
42 6ZL 


•> 


\ 


7^ 


L5 33L 
46 65L 




I 


75 


L5 5F 








00 IF 




Form l/l6 e x K^x) 


76 


40 F 




for small arg. 




22 57L 




77 


L5 199L 




\ 




LO 180L 




\ 


78 


40 8F 
50 176L 




■ 


79 


75 3F 




\ 




00 4f 




: 


80 


40 F 




: 




L5 8F 






81 


10 5F 




! 




LO F 




1 


82 


40 8F 
50 178L 






83 


75 8F 








00 4F 


4 


/ 


84 


40 8F 
26 87L 






85 . 


L5 86L 
46 9^L 






86 


50 138L 
26 88L 




- 


87 


L5 42L 
46 9^L 




■« 



[ LOCATION 


ORDER 


nam pass 6 v 11 


88 


L5 9^L 
42 62L 








y ».; 


L5 7F 
40 9F 




Place l/l6 e~ x I Q (x) 
in 9F 




90 


L5 91L 
42 72L 


1 


i 




91 


26 5^L 










50 92L 




1/16 e~* I-^x) for 




92 


50 5F 
75 7F 




saall arg. 




93 


00 IF 
40 7F 


* 


1 




9* 


26 F 








95 


50 163L 
L5 L 
42 62L 


; 


•\ 




96 


L5 44L 
46 65L 








97 


L5 5F 
00 IF 


1 






98 

99 

: 

: 

100 


40 F 
22 ^7L 

L5 199L 
%4 169L 
40 10F 

50 176L 


; 






101 


75 3F 
00 3F 








102 


40 F 










50 F 




l/l6x-e x K,(x) for 




i 103 


75 IF 




► 






00 3F 


. 


small argo 




104 


40 F 
L5 10F 








105 


10 3F 










L4 F 









LOCATION 



ORDER 



106 

107 
108 
109 

110 
112 
112 
113 

115 
116 

117 
118 

119 
120 
121 
122 
123 



kO F 
50 F 
75 178L 
00 5F 
kO 10F 
22 F 
L5 6F 
k2 111L 
50 136L 
50 191L 
L5 190L 
10 F 
66 5F 
00 2F 
S5 132L 
kO 3F 
7J 203L 
50 114L 
26 (Rl) 
i4-0 4F 
L5 H9L 
42 121L 
L5 113L 
k6 131L 
4l 199L 
L5 161L 
40 162L 
50 219L 
L5 3F 
40 2F 
50 2F 
7 J F 
IA 199L 
kO 199L 
F5 162L 
kO 162L 
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J 



► Form 3.75/X 



/ 



Form l/V^" 



Terms for 



large arg. 



LOCATION 


ORDER 
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12k 


52 127L 
50 5F 






125 


7 J 2F 
kO 2F 






126 


F5 121L 








1*2 121L 






127 


26 121L 
L5 199L 


4 


J 


128 


IA 218L 




\ 




kO F 




for large arg. 


129 


50 F 
75 ^F 




1 


150 


10 kF 
kO 7F 




I l/l6 e" x I Q (x) 


151 


26 F 
50 22 7L 




J 1/16 e" x I^x) 


152 


L5 151L 
k2 121L 




' 


155 


L5 7F 








ho 9F 




/ 


13k 


L5 110L 
k6 151L 






155 


26 118L 
50 206L 






156 


50 3F 
75 204L 






137 


^0 3F 
26 1^5L 




2/X 


158 


50 F 
50 191L 






139 


L5 I69L 
66 5F 






lto 


S5 150L 








40 5F 




2/X 


liH 


10 IF 






■ 


50 iklL 

















LOCATION 


ORDER 
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ite 


26 (Rl) 
40 4F 




i/ix 


14} 


16 135L 
42 62L 


■ * 


V 


144 


L5 1^9L 
42 153L 






145 


L5 140L 
46 65L 






146 


L5 3F 






; 


40 IF 






147 
148 


40 2F 

41 199L 
L5 200L 




l/l6 e X K^x) 


: 

149 


40 201L 
26 62L 




► 1/16 e x KqCx) 


• 


50 154L 






150 


L5 199L 
L4 205L 




for large arg. 


151 


40 F 








50 F 






152 


75 4f 
10 3F 






153 


40 10F 
26 F 






154 


L5 16OL 
42 62L 






155 


L5 10F 
40 8F 






156 


L5 157L 
42 153L 






157 


26 145L 
50 108L 






158 


L5 23% 
40 8F 




1 


159 


L5 256L 








40 10F 


■ > 

4 4- 


; 



LOCATION 


ORDER 
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160 


22 108L 




Return to master 




50 212L 




waste 


161 


LL 4095F 








IX kO&QF 




-8 


162 


00 F 
00 F 






I63 


00 F 




s 




00 1098 632^ 2500 J 




■164 


OOF 








00 6^37 3586 250 J 


Coeff. for 


165 


00 F 




» 




00 1885 6167 


500 J 


iito 


166 


00 F 




small arg« 




00 3323 ^162 250 J 




167 


00 F 

00 3769 15000 J 




168 


00 F 






00 ^051 3750 J 


1 


169 


00 F 








00 5000 0000' 0000 J 




170 


00 F 








00 7721 5720* 000 J 




171 


80 F 1 






00 6636 0710 5000 J 


:oeff. for 


172 


80 F 








00 9092 1551' 0000 J 




173 


80 F 








00 990^ 0299 4 0000 J 


Small arg. 


17^ 


80 F j 






00 999 1 * 4798 0000 J 




' 175 


80 F 






00 9999 7657 0000 J 




176 


00 F 








00 F 


_J 





LOCATION 
177 

178 

179 
180 
181 
182 

183 

181* 
185 
186 
187 
188 
189 
190 

191 
192 

193 



ORDER 

00 F 
00 F 

00 F 

00 F 

00 F 

00 1143 5526 312 J 

00 F I 

00 5772 1566 0000 J 

00 F 

00 4227 8420 0000 J 

00 F 

00 2306 9756 0000 J 

00 F 

00 3488 5900 000 J 

00 F 

00 2626 980 000 J 

00 F 

00 1075 00000 J 

00 F 

00 7400 000 J 

00 F 

00 3465 7359 0280 J 

00 F 

00 F 

00 F 

00 F 

00 F 

00 9375 0000 0000 J 

00 F 

00 F 

00 F 

00 6250 0000 000 J 

00 F 

00 2197 2643 1200 J 



HOTSS 



PAGE 11 



V U 



1 



Coeff. for 
KqCx) 

small arg. 



L LOCATION 


ORDER 


NOTES 
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194 


00 F 










00 1931 2140 0000 J 


Coeff 




195 


00 F 1 








00 ?5fe2 1825 000 J 


Uorl(x) 




196 


00 F 


\J 


- 


! 


00 1662 3325 000 J 


small arg. 




197 


00 F 1 

00 225^ 8000 00 J 






198 


00 F 1 








00 2863 3125 J 


/ 




199 


00 F 
00 F 








200 


LL 4095F 
LL 4090F 




-0 




201 


00 F 
00 F 








202 


00 F 










00 13^0 3565 603 J 






203 


00 F 1 

00 2666 6666 6667 J 






204 


00 F 1 

00 5333 3333 3334 J 






205 


00 F 1 

00 6266 5707 0000 J 






206 


80F 










00 9608 3* 


21 0000 J 


Coeff. for 




207 


00 F 










00 1094 7£ 
■ 


40 000 J 


)yx) 




208 


80 F 






1 


00 9946 8777 0000 J 
■ 


large arg. 




209 


00 F 






i 


00 2939 3600 00 J 






210 


80 F 1 






j 


00 9987 4230 0000 J 







LOCATION 
211 



212 
213 
21^ 

215 
216 

21? 
218 
219 
220 
221 
222 
223 
22^ 



226 
227 



ORDER 

00 F 
00 2660 
00 F J 

00 117^ 

80 F I 

00 9817 

00 F I 
1 

00 7521 

80 F I 

00 9960 

00 F 

00 1628 

80 F 

00 9996 

00 F I 

00 3989 

00 F I 

00 1328 

00 F I 

00 2253 

80 F I 

00 998^ 

00 F 1 

00 9162 

80 F I 
00 979^ 
00 F I 
00 2635 

80 F 

00 9835 

00 F j 

00 3923 

80 F 



00 960: 



1+000 J 
9309 3000 J 
2190 0000 J 
3^00 00 J 

9823 5000 J 
0700 00 J 
5877 5000 J 

1*228 0000 J 
5917 000 J 
1870 00 J 
21+35 1000 J 
8080 00 J 
2293 7000 J 
5372 000 J 
2367 1000 J 
7670 00 J 
1975 8000 J I 



[J 



NOTES 
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^^fe 



Coeff . for 



► K^x) 



large arg. 



Coeff. for 



^ 



large arg* 



LOCATION 
228 

229 

230 

231 

232 

233 

235 
236 

237 
238 



ORDER 
80 F 

00 9963 7981 7000 J 

00 F J 

00 1638 OlkO 00 J 
80 F I 
00 9896 8kk5 0000 J 
00 F J 

00 2282 9673 000 J 

80 F J 

00 9710 k6&J 9000 J 
00 F J 

00 1787 6535 000 J 

80 F J 

00 9957 99^1 3000 J 
00 F j 

00 5259 8013 292 J 

00 F I 

00 6^59 2302 791 J 

L3 188L I 

36 9L 

22 3L 

50 F 

(Rl) OOK 

(S3) OOK 

CR5) ook 

(Sk) OOK 



EOTES 



PAGE Ik 



V 11 



Coeffo for 

fl x (x) 
large argo 



239-247 
2^8-261 
262»278 
279=299 



